Efficiency of energy transfer in a light-harvesting system under quantum coherence 
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We investigate the role of quantum coherence in the efficiency of excitation transfer in a ring-hub 
arrangement of interacting two-level systems, mimicking a light-harvesting antenna connected to a 
reaction center as it is found in natural photosynthetic systems. By using a quantum jump approach, 
we demonstrate that in the presence of quantum coherent energy transfer and energetic disorder, the 
efficiency of excitation transfer from the antenna to the reaction center depends intimately on the 
quantum superposition properties of the initial state. In particular, we find that efficiency is sensitive 
to symmetric and asymmetric superposition of states in the basis of localized excitations, indicating 
that initial state properties can be used as a efficiency control parameter at low temperatures. 



I. INTRODUCTION 

Solar energy conversion in photosynthetic bacteria 
relies on sophisticated light-harvesting (LH) antennae 
which capture photons and then transfer the electronic 
excitation to a molecular complex which serves as a re- 
action centre (RC). There charge separation takes place 
and chemical energy storage is initiated (e.g., see Ref. 
1). Some of the harvesting complexes (LHl) and the 
RC are closely associated and form a core unit to ensure 
an efficient pathway for the transfer of excitations com- 
ing from peripheral antennae (LH2). This transfer takes 
only a few hundred picoseconds and is performed with ex- 
traordinarily high efficiency: most of the absorbed pho- 
tons give rise to a charge separation evenii. The precise 
mechanisms underlying such high efficiency remain elu- 
sive despite numerous studies on the subject. In partic- 
ular, whether quantum coherence plays any role on pro- 
moting the efficiency is still ambiguous. Some works in- 
dicate that it will induce higher excitation-transfer rates^^ 
while others argue that this may not necessarily be the 
case^. Remarkably, recent experimental and theoreti- 
cal works^iS indicate that long-lasting electronic coher- 
ence can indeed influence the excitation transfer dynam- 
ics in photosynthetic complexes. For instance, quantum 
beats associated with electronic coherence in the Fenna- 
Matthews-Olson (FMO) complex of green sulfur bacteria, 
which connects a large LH to the RC has been reported 
by Engel et al.'^. Also, coherence among electronic states 
of closely associated pigments in the RC of purple bac- 
teria has been reported by Lee et al.^ . Furthermore, the 
excitation transfer in organic dendrimers has attracted 
significant attention through the prospect of creating ar- 
tificial photosynthetic systems^. One of the key obser- 
vations in these artificial systems is the evidence of co- 
herent energy transfer mechanisms. These experiments 
therefore open up the possibility of exploring in detail the 
interplay between quantum coherence and the efficiency 
of natural and artificial LH systems. 

In this work we consider as model system a ring-hub 
arrangement of interacting two-level systems, represent- 
ing a LHl-RC core unit as in purple bacteri a^'^'^" , and 
use the quantum jump approac h^^i^^ to provide a simple 
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FIG. 1: Schematic of the LHl-RC core of purple bacteria 
Rodobacter Sphaeroides. (a) Arrangement of the 32 Bacte- 
rioclorophils (BChl) surrounding the RC. The RC has two 
accessory BChl and two acceptors forming a special pair re- 
sponsible for charge separation, (b) Diagram of the induced 
dipole moments in (a). The arrows indicate the dipole mo- 
ment directions corresponding to data taken from Xhu et al^. 
(c) Toy model: the RC is assumed to be a single two-level sys- 
tem. 



picture of how the quantum superposition properties of 
the initial state of the excitation relate to the efficiency of 
transfer from the LHl complex to the RC when coherent 
energy transfer and static disorder of single-site energies 



dominate. The quantum jump approach 



proves to 



be particularly suitable to describe excitation dynamics 
in this situation because the density matrix elements can 
be calculated exactly from the no-jump evolution, given 
that there is only one excitation at most in the system 
and that only local dissipation rates (or charge separa- 
tion rates) are assumed^. Considering initial states as 
superpositions of single-site excitations, we find that the 
efficiency profile depends both on the symmetry proper- 
ties of such superposition states, and on the number of 
sites among which the excitation is initially delocalized. 
In particular, our results show a non-trivial interplay be- 
tween excitation delocalization and efficiency as there can 
be an optimal delocalization length for which efficiency 
of transfer reaches a maximum and transfer time a min- 
imum. Such behaviour is robust to the presence of en- 
ergetic disorder. The plan of this paper is as follows. 
The next section outlines the model of energy transfer at 
very low temperatures as well as how we use the quan- 
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turn jump approach to calculate the main characteris- 
tics of our photosynthetic system. Section III discusses 
our results both for a toy model and for a more detailed 
Hamiltonian describing the LHl-RC complex. 



II. COHERENT EXCITATION TRANSFER 

We consider a system of M donor pigments surround- 
ing a RC with N acceptors (sec Fig. [T|) described by the 
Hamiltonian H — Hq + Hj. Labeling the donors from 1 
to M, and the acceptors from M+1 to M + N, the single 

particle Hamiltonian is Hq = X^^l^i^ ^a<^a'^a where Eq, 
is the excitation energy of pigment a, and the interaction 
Hamiltonian reads: 

M M+N M M+N 

Hi ^ 7jcV,c 4- ^ JjkVjk+ ^ QcrVcr ■ 

j — 1 c — A/ + 1 J— l;fc>j c=M + l;r>c 

Here Vat ^ <^t<^b + '^t'^a where (t+(^) is the Pauli op- 
erator for a two-level system, and 7jc, Jjk and gcr are 
the donor-acceptor, donor-donor, and acceptor-acceptor 
couplings, respectively. We are interested in the low tem- 
perature regime where static disorder dominates and dy- 
namical effects can be neglectedi^. In order to account 
for static disorder we treat as having a random com- 
ponent: Ca = Ea{0) + whcrc Ea{0) is the ensem- 
ble average value, and 5Ea is the energy disorder at site 
a given by a Gaussian distribution with zero mean and 
standard deviation a. We assume that the open system 
dynamics is dominated by two incoherent processes: the 
excitation can be dissipated in a donor or it can induce 
charge separation at a site in the RC. Such dynamics of 
our coherent photosynthetic core can be described by the 
Lindblad master equation (h — 1): 

^p = -i[//,p] + i^[2yl<,p^I-AU<,p-p^U<.] , (2) 

where the commutator generates the coherent part of the 
evolution and the action of each operator = \/2Ta<J^ 
accounts for a "jump" process associated either to dis- 
sipation of excitation in a donor i.e = F with a — 
or for a charge separation event at an accep- 
tor of the RC i.e = k with a = M + 1, . . . M + N . 
Here we have assumed identical dissipation rates for the 
donors and identical charge separation rates for the ac- 
ceptors at the RC. This formalism can be extended to 
include other incoherent processes. 

In order to solve equation ([2]) we follow the quantum 
jump approacbii and re- write equation ^ as 

p = -t{H,ondP - pHl^d) + "^"P^^o. (3) 

a 

with 

M N 

Hcond = - «F ^ (T^aJ - in ^ cr+o-;: , (4) 

j=l c=M+l 



In this description, the excitation dynamics can be in- 
terpreted in terms of quantum trajectories where the 
system follows a no-jump evolution associated to the 
non-hermitian Hamiltonian Hcond, interrupted by a sin- 
gle stochastic collapse of the system to its ground state 
in the event of either dissipation or charge separation, 
with probability Pa = iT:[AapA\^. The no-jump tra- 
jectory conditioned on no-decay-occur is described by 

Pcond{t) = -i{HcondP - pHl^nd)- particular, if the 
initial state is pure, i.e. |^(0)), the state remains pure 
(but unnormalized, i.e. dissipative) in the no-jump tra- 
jectory and becomes |*corirf(0) = exp{-iHcondt)\'if{0)). 

We now demonstrate that given that a single excita- 
tion is present in our photosynthetic core, the dynam- 
ics of all the density matrix elements can be calculated 
exactly knowing only Pcond{t). Notice that Hcond pre- 
serves the number of excitations i.e. [i/cond, A/"] = with 

N = Y^^a=i' ^a'^a- Heucc, for a single excitation the 
density matrix dynamics is restricted to the subspace 
of single-excitation states plus the ground state. We 
choose the basis given by 5 = {|0), {|c)}}, where 

|0) = |0i . . . O^f ; Qm+1 ■ ■ - ^m+n) is the state with all the 
pigments in their ground state, and 

\j) = |0i . . . Ij . . .OAf;OM+i ■ • -OAf+jv) 

|c) ^ |0i . . .OAf;OAf+l . . . Ic • • •OAf+A') 

are states in which only the jih. donor (or cth accep- 
tor) is excited. The labels after the semicolon in each 
ket refer to the acceptors at the RC. Let us denote 
the density matrix elements pki{t) = {k\p{t)\l) for any 
pair of states |fc),|0 G S. Notice that the second term 
of Eq. Q satisfies {k\Aap{t)A'^c,\^) = for all states 
except when \k) = \l) — |0). Therefore, for single- 
excitation states pM{t) = {k\pcond{t)\l) , while pooit) — 
EjO\^Mt)AUO) = 2E„r„(a|p(i)|a). Now, since 
p{0) = pcondiO), then pooit) = 2J2a^a{a\pcondit)\a). 
This demonstrates that the dynamics of all density ma- 
trix elements can be entirely calculated with Pcond(t) and 
hence our claim. 



A. Efficiency and Transfer times 

With the above formalism we can now focus on the 
main features of our coherent LHl-RC core. As it is 
described in the review by Sener et alM-, of particular 
interest are the efhciency, which is given by the probabil- 
ity of an excitation to be used for charge separation as 
opposed to being dissipated, the average transfer time of 
excitation to get trapped by the RC, and the excitation 
lifetime after the initial absorption of a photon. 

Let us denote the initial state '^q. Clearly, the proba- 
bility that the excitation is still in the system at time i, i.e 
no-jump probability is P{t;'^o) = J2k'^^ {^'\Pcond{t)\k) 
while w{t; '^q) = poo{t) is the probability density that a 
'jump' (charge separation or dissipation) occurs between 
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[t, t + dt) and it reads 

M M+N 
W(t;*0) = 2r^{j\pcond{t)\j) +2k ^ {c\pcond{t)\c) 
i = l c-M+l 

= wnit;'i/o) + WRc{t;<bo) ■ (5) 

Here WD{t]'^o)dt is the probability that it is dissipated 
by any of the donors in [t, i + dt) while WB.cit]'^o)dt 
is the probability that the excitation is used for charge 
separation at the RC. Notice also that w(i;^'o) = 
-dP{t]^f))/dt leading to w{t;^!Q)dt = 1 which im- 
plies that the excitation will eventually either be dissi- 
pated or trapped in the RC. In particular, for pure initial 
states of the form ^-o = ''j(O)li) +EcSf+i &c(0)|c) 

with Y.f=i hm'' + Ecilvf+i l^c(0)P = 1, we have that 
Pcond{t) = |*cond(i))(*cond(i)| with the unnormalized 
conditional state given by 

M M+N 
\^cond{t)) ^J2^,{t)\j) + ^M0|C). (6) 
j=l c=M+l 

The monotonically decreasing norm of this state gives 
the no-jump probability P{t;^o) =|| l^'cond) |P, while 
WD{t;^o) = 2Tj:f^,\b,it)\^, and 

M+N 
c=M+l 

We therefore define the efficiency (rj) of energy transfer 
to the RC as the total probability that the excitation 
is used in charge separation. The transfer time (tf) is 
the average waiting-time before a jump associated with 
charge-separation in the RC, given that the excitation 
was initially in the LHl ring. The excitation lifetime (r) 
is the average waiting-time before a jump of any kind 
occurs: 

T)= i dtwRc{t--io) , tf = - I dttwRcit--^o) , 
Jo "l Jo 

poo 

r = dttw{t; *o) • (8) 
Jo 



III. EFFICIENCY CONTROL MECHANISMS 

In the classical descriptionMii^, where incoherent ex- 
citation transfer is assumed, typical efficiencies of energy 
transfer tend to be near unit. This is due to a separa- 
tion of the dissipation (ns) and the excitation transfer 
and charge separation time scales (ps). We shall shortly 
show that using the same parameters for single-site ener- 
gies, electronic couplings, dissipation and charge separa- 
tion rates, the efficiencies obtained under coherent trans- 
fer are much lower and strongly dependent on the initial 
state of the excitation. We find that the initial rela- 
tive phases between localized excitation states \j), and 
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FIG. 2: Numerical results for (a) r/, (b) tf, and and (d) r 
versus m, for the toy model with three different interaction 
mechanisms. For nearest- neighbour interactions (O) the cou- 
pling is lOOmeV, while for the pairwise case it is lOmeV 
and equals the average dipole-dipole coupling (•). In each 
case, the donor- RC coupling 7 equals ImeV, F = 1 ns~^ and 
K = 4 ps~^. (c) No-jump probability for the case of dipole- 
dipole interactions, as a function of time and for different m 
values. 



the number of donors among which the excitation is ini- 
tially delocalized, can act act as efficiency control mech- 
anisms. In what follows, we first consider a simple model 
for which analytical solutions can be obtained and then 
we calculate efficiency and transfer times with the model 
Hamiltonian given by Hu et al^^. 



A. A toy model for the LHl-RC complex 

The simplest model for which analytical solutions can 
be obtained corresponds to the RC taken as a single two- 
level system i.e. TV = 1, on resonance with the M = 32 
donors in the LHl ring (see Fig. [IJc)). Later we shall 
show that the main qualitative behaviour observed in 
this situation also applies to a model featuring the de- 
tailed structure of the LHl-RC complex in purple bacte- 
ria. We consider initial states in which the excitation is 
delocalized among donors, i.e. ^I^o = ^jLi ^i(O)b) with 
X]j=i — 1- The unnormahzed state of Eq. ^ be- 

comes \^cond{t)) = Ejii^jWIj) + bM+i{t)\M + l) sat- 
isfying the equation dl^i condit)) / dt = -iHcond\^condit)) 
which leads to a set of first-order coupled differential 
equations for the complex amplitudes bj{t) and bM+i{t). 
In Appendix A we show that when the system's dy- 
namics is invariant with respect to exchange of donors 
in the antenna, one can find analytical solutions for 
bM+i{t). Let us define the effective interaction between 
a donor j and the rest of pigments in the LHl ring as 
Aj — J2k '^jk the sum of all the coupling strengths 
between donor j and any other pigment in the LH ring. 
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The system's dynamics is invariant with respect to donor 
exchange when both all donor-RC couplings are iden- 
tical, i.e. 7jc = 7, and the effective interaction be- 
tween a donor and the rest of pigments in the ring 
are also identical for all donors, i.e. Aj = A for all 
j. Under these conditions bM+i{t) satisfies the differen- 
tial equation bM+i{t) + XbM+i{t) + YbM+i{t) = with 
X = {k + r + iA) and Y ^ iMj^ + K(r + iA) (see 
details in Appendix A). For the initial condition where 
&M-i-i(0) = i.e. excitation is initially in the ring, we 
find 

M 

\bM+i{t)\' = m\Y.^Ao)\' ■ (9) 

Here T{t) = 472e-(^+'")*|sin(f7i/2)| Vjr^P and ft is the 
complex frequency that determines the timescale of co- 
herent oscillations i.e. fl = y/4AI^^ — (F — k + iA)^. 
Note that A identical for all donors does not imply that 
the pair couplings Jjk need to be identical for all possi- 
ble pairs. Hence, analytical solutions can be found for 
three different mechanisms of interaction between the 
donors: (i) nearest neighbours with Jjj+i — J/2 (ii) 
pairwise interaction with Jjk = J for all (j, k) pairs and 
(iii) dipole-dipole interactions of the form Jjk = J/f^k 
where ijk is the relative position vector between the in- 
duced dipole moments of donors j and k. From equation 
([7|) the probability density of having charge separation 
becomes w_Rc(i;*o) = 2k|6a/+i(0P; then we obtain an 
expression for the corresponding efficiency: 

/•oo 

T] = 2k\Bo\^ / T{t)dt (10) 

with Bq = Equation (fTn|) is the main re- 

sult of this paper. It shows that the efficiency of transfer 
depends on the quantum coherence properties of the ini- 
tial state as it becomes proportional to |-BoP i-e. the 
amplitude of probability, and not just the probability, 
that the excitation is initially in the LHl. Therefore, 
the efficiency profile is sensitive to symmetric and asym- 
metric superpositions of localized excitation sates \j) i.e. 
it depends on the initial relative phases between states 
From equation pl}|) one can conclude that symmet- 
ric delocalized excitation states yield an increase in 77, 
while some asymmetric states could be used to limit or 
even prevent the transfer, i.e. 77 = 0. Unless otherwise 
stated, we henceforth consider symmetric initial states 
of the form = (l/\/"i) X^jli b) where m < Af is 

the number of donors among which the excitation is ini- 
tially delocalized. We denote m the delocalization length. 
For these symmetric states |i?oP — ™ and hence 77 oc m 
as shown in Fig. ^a). This figure also shows that the 
efficiency gradient depends on strength of the interac- 
tion between one donor and the rest, which is quanti- 
fied by A. We have chosen 7 to be the same for all 
these situations, but J has been taken to be such that 

Anearest < Adipole — Apairwise- For a fixcd 771, Tj rCachcS 



higher values in the case of nearest-neighbor couplings, 
while it achieves similar values for dipole-dipole and pair- 
wise interactions. According to these results, interaction 
among donors limits the efficiency: the stronger the effec- 
tive interaction between one donor and the ring, the lower 
the efficiency will be. This phenomenon seems to resem- 
ble the 'entanglement sharing' dynamics in the context 
of a central spin coupled to a spin-bath~.. In our case the 
LHl complex can be seen as a spin bath for the RC. Daw- 
son et ali^ have discussed that interaction between bath 
spins translates to entanglement among them, and since 
entanglement cannot be shared arbitrarily among several 
particles, interaction among spins in the bath limit en- 
tanglement between the central spin and the bath, and 
therefore may limit the efficiency of transfer. A discus- 
sion of efficiency in terms of entanglement is beyond the 
scope of this paper and the work in this direction will 
be presented elsewhere. Interestingly, tf in this simple 
model turns out to be independent of 771, as can be de- 
duced from Eq. (7). Therefore, for the symmetric initial 
states considered, tf depends mainly on the mechanism of 
interaction, as can be seen in Fig. ^h). The decay-rate 
of P{t; l^'^j)) increases with 777, as shown in Fig. [DJc). 
Correspondingly, the excitation lifetime r decreases as 
shown in Fig. [2td). The three situations satisfy tf < t, 
where the equality holds for the initial state in which the 
excitation is symmetrically delocalized among all donors. 



B. A detailed model for the LHl-RC complex. 

We now apply the above formalism to the effective 
Hamiltonian for the LHl-RC interaction given in refer- 
ences 9 and lO^iiiS. In this case the RC has a special pair 
of BChl responsible for the charge separation, i.e. N = 2 
acceptors, and two more accessory BChl molecules which 
do not participate in the charge separation processi^ (see 
Fig. [1]) . The effective Hamiltonian is of the form given in 
equation ([T]) but with certain particularities. First, the 
pigments at the RC are off-resonance with the donors. 
Second, the interactions between adjacent molecules are 
quantified by two different constants i.e. J2j.2j+i = 1^1 
and J2j,2j-i — ^2 which are derived through quantum 
chemical calculations^' Third, the coupling between 
non-neighbouring donors corresponds to a dipole-dipole 

interaction of the form J^k = - ^^^^'■•'^^i^^^"'^"^ 

where jlj is the transition dipole moment of the j*'* donor 
and fjk is the relative position vector between donors j 
and k. The directions of /ij have been taken from Hu et 
al.^'^, and a top view of the dipole representation of the 
LHl-RC core is shown in figure[ljb). For the initial con- 
dition of the excitation in the LHl complex we consider 
both symmetric initial states l^'^) = (l/\/"i) X^jli b) 
with total amplitude of probability l-BoP ~ and asym- 
metric states l^-^") = (l/V^)EJLi(-l)-'b') satisfying 
I-SqP = 0. Since the system's dynamics is not invariant 
with respect to pigment exchange, for each delocalization 
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FIG. 3: Numerical results for initial symmetric states of 
the form |*^) = (l/Vm) b) satisfying \Bo\^ = m. 

(a) Efficiency and (b) transfer time versus m for different 
cr— values. Results shown are averaged over an ensemble 
of 1000 aggregates for each standard deviation a. Single- 
site energies and electronic couplings for the LHl-RC core 
have been taken referencei^. For donors in the LHl com- 
plex -Ej(O) — 12911 cm~^ while for the special pair at 
the RC Ec{Q) = 12748 cm"\ and for the accessory BChl. 
Ea{Q) = 12338 cm"^ In all cases F = 1 ns"^ and « = 4 ps"^ 



FIG. 4: Numerical results for initial asymmetric states of the 
form = (l/V^)Er=i(-l)'b') satisfying |Bo|' = 0. 

(a) Efficiency and (b) transfer time versus m for different 
a— values. Results shown are averaged over an ensemble 
of 1000 aggregates for each standard deviation o. Single- 
site energies and electronic couplings for the LHl-RC core 
have been taken from reference^S. For donors in the LHl 
complex Ej{Q) = 
the RC Ec{Q) = 



12911 cm while for the special pair at 
12748 cm~^, and for the accessory BChl. 



Ea{Q) = 12338 cm~^ In aU cases F = 1 ns"^ and k = 4 ps" 



length m, we calculate the average efficiency rjm = {i]) 
where the average is taken over all possible states for 
which the excitation is delocalized among m consecutive 
donors. Also, when energetic disorder is considered, for 
each value of standard deviation cr, the efficiency corre- 
sponds to the ensemble average over 1000 realizations of 
disordered H. We denote this average efficiency {rim) a 
and the corresponding average transfer time {tf)a- An 
estimate of cr = 30cm~^ for the standard deviation of the 
diagonal disorder distribution in the LHl has been pro- 
vided in the literature^^. Therefore we have carried out 
calculations for cr— values ranging from to ISOcm"-'^ and 
the results for symmetric and asymmetric initial states 
are shown in figures [3] and [J] respectively. For symmetric 
states and in the absence of energetic disorder i.e. cr = 0, 
the behaviour of the efficiency and transfer time are very 
similar to that in the toy model: (7?m)o-=o increases lin- 
early with m (see Fig. [31[a)) while {tf)a=o decreases (see 
Fig. Efb)), achieving maximum efficiency and minimum 
transfer time for the fully delocalized situation. For dis- 
order distributions corresponding to small values of cr 
i.e cr < SOcm""'^, there are symmetric states which ex- 
hibit efficiency improved in comparison to the situation 
where no disorder is considered. Such states correspond 
to those for which m < 10 as it can be seen in figure 
[3]Ja). However, the efficiency values are clearly decreased 
with increasing cr. Conversely for the asymmetric states 
satisfying |i?oP = 0, the efficiency is a non-monotonic 
function of m, indicating that there is an optimal delo- 
calization length for which {rjm) a has a maximum and 
for which {tf)a- has a minimum as it is shown in figure 



m The optimal delocalization length is around to = 10 
as it can been see in figure Ufa). Such behaviour is ro- 
bust to the presence of energetic disorder, and indeed im- 
proved as the efficiency values are larger with increasing 
cr. It is also worth noting that in both cases, symmetric 
and asymmetric states, the efficiencies obtained are lower 
than those given by a classical calculation with rate equa- 
tions derived for the same single-site energies, electronic 
couplings, and dissipation and charge separation rates 
here considered^. 

The above results suggest that efficiency can therefore 
be used as an indicator for coherent energy transfer. In 
particular, two-dimensional spectroscopy techniques re- 
cently developed to study coherence dynamics in pho- 
tosynthetic systems*'^ may be used to create and probe 
quantum superposition initial states. For instance, an 
optically allowed state of LHl is the completely delocal- 
ized asymmetric stated. For very low values of energetic 
disorder and at very low temperatures, the efficiency of 
transfer from such state is nearly zero. At a slightly 
higher temperature, the excitation will become less de- 
localized and as such the efficiency would increase (c.f 
Fig. [¥]). In other words, an experiment measuring the ef- 
ficiency of an LHl-RC core under various temperatures, 
but still within a low temperature regime, could serve to 
ascertain the extend of coherence in energy transfer. In- 
terestingly, some experimental works have indicated that 
in thermalized LH2 complexes the excitation may be co- 
herently delocalized over just a few donors of the B850 
ringi^ while at very low temperatures it can be fully de- 
localized over the whole ring^". Unfortunately, no such 
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investigation has been reported on the LHl. 

In conclusion, we have presented a formalism to study 
the role of quantum coherence in photosynthetic units 
that exhibit coherent energy transfer mechanism. Our 
results open up experimental possibilities to investigate 
and exploit such coherent phenomena in artificial and 
natural systems harvesting light. 
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with 



M 



M M 



git) = ^ h, {t)-iY.T. J^^^^ (^) 



j=i k^j 



Defining the effective coupling between one donor and the 
LHl ring as Aj = Sfc^j '^jki assuming it identical 
for all the donors in the LHl i.e. Aj = A, we have 
g{t) = 7(ir - A) J2jLi bjit). Now from equation 

we get Y.f=i bj{t) = i7"^(^A/+i(t) + nbM+iit)) and thus 
we arrive to 



APPENDIX A: DETAILS OF THE ANALYTICAL 
SOLUTIONS FOR THE TOY MODEL 

As described in section HI (A), in the simplest situ- 
ation where the RC is taken as a single two-level sys- 
tem on resonance with M donors, we are able to find 
analytical solutions for the probability density of hav- 
ing a charge separation event at the RC, and there- 
fore we find an analytical expression for the efficiency 
of transfer in this toy model (see Eq. (fTO|l ). Here we 
give the details of this calculation. The unnormalized 
state \'^cond{t)) given in equation ^ satisfies the rela- 
tion d\'i>cond{t))/dt = -iHcond\^cond{t)) which leads to 
the following set of first-order coupled differential equa- 
tions 



&M+i(t) -t'^J,kbk{t) - Tb,(t) (Al) 



M 



6jv/+i(i) = ~^y^7jfej(^) - nbM+i(t) 



(A2) 



Assuming that all donor-RC coupling are identical i.e. 
7j = 7 for all j = 1 . . . 32, we have 

bM+lit) + KbM+lit) + M-1%M+l{t) - Q{t) 



bu+i (t) + XbM+i[t) + YbM+i (i) = , (A3) 
with A = (k r zA) and y = 4M72 + k{Y + iA). 
The solutions of the above differential equation are of 
the form 



M 



bM+i(t) = f{t)Y^b,{0)+g{t)bM+m (A4) 



with 



fit) = _2z7e-^*/^^^^^^ 



9{t) = e 



(F - K + iA) 



sin(m/2) +cos(m/2) 



and — y/4Mj'^ — (F — k -I- zA)^. The above solution is 
invariant with respect to exchange of any pair of two-level 
systems. It is also worth noting that the above formal- 
ism allows to find a closed expression for the collective 
amplitude of probability Bit) = J2 
further simplifications of this system 



^j ' (^) which leads to 
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